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Abstract. - We present a theoretical and experimental analysis of the transverse monopole 
mode of a fast rotating Bose-Einstein condensate. The condensate's rotation frequency is 
similar to the trapping frequency and the effective confinement is only ensured by a weak quartic 
potential. We show that the non-harmonic character of the potential has a clear influence on 
the mode frequency, thus making the monopole mode a precise tool for the investigation of the 
fast rotation regime. 



The investigation of rotating gases or liquids is a central issue in the study of superflu- 
ids |1I2| . During the recent years, several experiments using rotating atomic Bose-Einstein con- 
densates have provided a spectacular illustration of the notion of quantized vortices |3I4I5I6| . 
Depending on the rotation frequency of the gas, a single vortex or a regular array of vor- 
tices can be observed experimentally. When the rotation frequency is increased to a very 
large value, a new class of phenomena is predicted, in connection with quantum Hall physics 
|7ll51l^l lTUlllllll2iri3lll4| . For a gas confined in a harmonic potential, the fast rotation domain 
corresponds to stirring frequencies Q of the order of the trapping frequency u>± in the plane 
perpendicular to the rotation axis (hereafter denoted z). From a classical point of view, the 
transverse trapping and centrifugal forces compensate each other for this stirring frequency, 
and the motion of the particles in the xy plane is only driven by Coriolis and interatomic 
forces. This situation is similar to that of an electron gas in a magnetic field, since Lorentz 
and Coriolis forces have the same mathematical structure. 

In order to approach the regime of fast rotation two paths are currently being explored. 
The first approach is implemented in a pure harmonic potential and is based on evaporative 
spin- up, i.e. the selective removal of particles with low angular momentum |15j . The stirring 
frequency can then be raised close to u>± (O = 0.993 u>± was reached in ^H]). The second 
approach, which is followed here, consists in adding to the quadratic confinement a small 
positive quartic potential, which ensures that the particles will remain confined even when £1 
exceeds uj± |17I18I19I20I21I22] . We have recently proven that this method can be successfully 
implemented and we have mechanically stirred a rubidium Bose-Einstein condensate up to 
n ~ 1.05 wj. |25] . 

All experiments performed so far (including the present one) are still deeply within the 
mean field regime, characterized by a number of vortices N v well below the number of particles 
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N. Nevertheless in order to prepare for future investigations of possible quantum-Hall-like 
states, one must design proper tools of investigation, such as the study of the eigenmodes 
of the rotating gas. Possible examples are the transverse quadrupole modes, which allow 
for a measurement of the z component of the angular momentum |23| . and the Tkachenko 
oscillations of the vortex lattice, recently observed in JHj- Here we report on the experimental 
and theoretical study of the lowest transverse monopole mode of an ultra-cold gas of rubidium 
atoms in the fast rotation regime. 

In a trap with axial symmetry along the z-axis, excitations can be characterized by their 
angular momentum along z (quantum number m z ) |24j . Here we are interested in excitations 
carrying no angular momentum (m z =0). In the limit lu z <C lo±_ which is of interest here, 
it is possible to identify m z = "transverse monopole" modes, which mostly affect the atom 
distribution in the xy plane. For a non-rotating condensate the lowest frequency for these 
m = transverse modes is oj mp ~ 2u>± (transverse breathing mode) [251126) . In the present 
work we show experimentally that the relation w mp ~ 2 to± remains valid for large rotation 
frequencies (as predicted by |27p. as long as the quartic term in the confinement is not 
significant. We also explore the region ~ where the quartic term plays an essential role. 
We compare the measured monopole frequency with that derived from a simple hydrodynamic 
model of an infinite, cylindrical condensate. Finally we discuss the structure of the mode, 
which exhibits a very particular behavior in the range f2 > lu±. 

Our 87 Rb Bose-Einstein condensate contains 3 x 10 5 atoms. It is produced by radio- 
frequency evaporation in a Ioffe-Pritchard magnetic trap, to which we superimpose a far-blue 
detuned laser beam. The magnetic trap provides a harmonic confinement with cylindrical 
symmetry along the z axis, with lu z /2t: = 11.0 Hz and uj^ /2n = 75.5 Hz. The laser beam 
adds a negative quadratic and a positive quartic potential in the xy plane. The quadratic 
term decreases the trap frequency uj±_ by ~ 15%. The quartic term allows to explore rotation 
frequencies around and slightly above the trapping frequency co±. It reads fcr 4 /4, with r 2 = 
x 2 + y 2 and k — 2.6(3) x 10~ n J/m 4 . The oscillation frequency of the condensate center- 
of-mass (dipolc motion) in the xy plane is Udp/^n — 65.6 (±0.3) Hz. Because of the quartic 
component of the trapping potential, cjdp is slightly larger than even for an arbitrarily 
small amplitude of the dipole motion. The condensate has indeed a finite radius R c — 6.5 fxxsx 
and therefore explores regions of space where the contribution of the quartic term is significant. 
Using a perturbative treatment of the quartic term we infer that u>d P — w± — 2k R 2 / (7mw 2 _), 
so that w±/(2tt) = 64.8 (±0.3) Hz. 

The procedure for setting the gas in rotation at a frequency f2 ~ to±_ has been described 
in detail in |23| . We apply an additional laser beam creating a rotating anisotropy in the xy 
plane for two successive phases, with respective durations of 300 ms and 600 ms. During the 
first phase the rotation at frequency Q ~ u)±/y2 resonantly excites the transverse quadrupole 
mode, and vortices subsequently penetrate the condensate. After a 400 ms relaxation period 
we apply the second stirring phase at the final desired frequency. This phase is followed by 
a 500 ms relaxation period. By measuring the transverse size of the condensate, we have 
checked that it is then rotating at a frequency close to the stirrer frequency O (within 2%) for 
the regime of interest fi < 1.05 u>±. In particular, for £1 < lj± we observe large vortex arrays 
involving up to ~ 50 vortices. During the stirring and relaxation periods, we apply radio 
frequency (r.f.) evaporative cooling. The r.f. is set 24 kHz above the value which removes 
all particles from the trap. Assuming f2 = ci>j_, so that the transverse confinement is purely 
quartic, this gives a trap depth Uo ~ 40 nK, hence a temperature T < 10 nK. 

We then excite the transverse monopole mode(s) of the gas by changing for a period r the 
intensity of the laser creating the quartic potential from / to I'. We choose r = 2 ms, which 
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Fig. 1 - Series of images of the condensate for various waiting times r after excitation. Row a) 
corresponds to a rotation frequency Q/(2tt) = 62 Hz, row b) to Q/(2n) = 64 Hz and row c) to 
n/(2rr) = 66 Hz. 



is short compared to the oscillation period of the monopole oscillation 27r/w mp = 7.5 ms. 
We then wait for an adjustable duration t before performing a time-of-flight expansion and 
absorption imaging. The images are taken along the rotation axis z so that we have access to 
the column density in the xy plane. 

Typical results are shown in figuren]for £l/2ir — 62, 64, 66 Hz. Each image corresponds to 
a destructive measurement of the atom density after a duration t = 25, . . . , 32 ms and a 18 ms 
time-of-flight. In order to check that the oscillation is in the linear regime, we explored values 
between and 4 for the ratio I' /I: we found that the frequency u)o is (within experimental 
uncertainties) independent of the excitation strength \I— I'\/I, and that the amplitude of the 
oscillation varies linearly with this strength. We have also varied the duration of the time-of- 
flight between 6 and 18 ms. We thus verify that all features reported here are independent of 
this duration, which corresponds to a mere scaling of the initial spatial density distribution. 

We have taken similar sequences of images for various stirring frequencies f2, with t varying 
from to 40 ms by steps of 1 ms. We extract from these images the transverse size R(t) of 
the gas. The variations of R(t) are well fitted by a single sinusoidal function, with frequency 
Wrap- The characteristic lifetime of the excitation is larger than 100 ms, so that its decay is 
negligible during the 40 ms period. The variations of function of are plotted in 

figure |2J An obvious feature of this plot is that ui mp varies very weakly for Q except in the 
vicinity of luj_. We now discuss in more details three relevant domains of this graph. 

(1) For a non rotating gas (f2 = 0) we find w mp /(27r) = 130.6 (±1.5) Hz, which is close to 
the well known result w mp ~ 2uj±_ = 2n x 129.6 (±0.6) Hz j2"5| . 

(2) When the gas rotates at a frequency fl not too close to uj±, the effect of the quartic 
potential remains small. In these conditions we find that the monopole frequency cuo stays 
approximately constant. At first sight this may seem a surprising result. For a rotating gas 
in a pure harmonic potential, the transverse confinement is reduced due to the centrifugal 
force. In the calculation of the equilibrium shape of the gas, the trapping frequency uj± is 
thus replaced by the weaker value {lo\ — J7 2 ) 1 / 2 . One might have expected that a similar 
replacement should be done also for the monopole frequency; this is clearly not the case. For 
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Fig. 2 - Transverse monopole frequency w mp (f2) normalized by cj mp (fi = 0) as a function of the 
stirring frequency Q.. For a non rotating condensate (Q = 0) we measure cj mp (0)/(27r) = 130.6 Hz. 
The vertical bars indicate the error given by the fitting routine. 

a cigar-shaped condensate in the hydrodynamic regime, it was proven in |27| that luq stays 
equal to 2cu± when f2 varies. In the case of an ideal gas described by classical mechanics, 
the equality o> mp = 2lu±_ also holds for any f2 as a consequence of the combined action of 
centrifugal and Coriolis forces. 

(3) When the rotation frequency Q approaches u>±, the contribution of the quartic term 
becomes more important and the monopole frequency o-> mp deviates significantly from 2oj j_. 
As shown in figure^ this deviation reaches ~ 8 % for f2/(27r) = 68 Hz, i.e. £l/ui± ~ 1.05. We 
have also plotted in figure |21 the prediction of a theoretical treatment of the monopole mode 
in the quartic + quadratic trap, neglecting the harmonic confinement along the z-dircction. 
The agreement between these predictions and the experimental data is satisfactory, given the 
difference in geometry between the experiment and the model. 

The starting point of our theoretical treatment is the set of equations for rotational hy- 
drodynamics [22], which we write in the laboratory frame: 

dtp = -V • (pv) (1) 
d t v = -V({U + gp)/m + v 2 /2) + v x (V x v) , (2) 

where p and v are the density and velocity field of the atom distribution. U = mu)\r 2 /2+kr 4 /4 
stands for the trapping potential and g — A^Y^ajra characterizes the strength of the atomic 
interactions (a is the scattering length). As explained in 29 2Tll3Tj] these equations are valid 
when one is interested in a phenomenon whose characteristic length scale is larger than the 
distance between vortices. 

The eigenmodes are obtained by linearizing 11121) around the rotating equilibrium solution 

v cq = ft x r gp cq = n - U - U ccn , (3) 

where /i is the chemical potential and U ccn = — m0 2 r 2 /2 the centrifugal potential. We elim- 
inate the velocity field to get a closed equation for the density variation dp, assuming an 
oscillation at frequency u>: 



-V ■ ( Pcq V[Sp}) = — {oj 2 - 4tt 2 ) [dp] 



(4) 
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In absence of a quartic term the eigenfrequencies for the m = modes are 



= 4f2 2 + 2n(n + l)(u)± - , ra positive integer, (5) 



and the corresponding eigenmodes are polynomials of degree n with respect to r 2 . For n = 1, 
we recover in particular u>i = uj mp = 2 uj± . 

When Q approaches u±, all transverse modes m = become degenerate. Such a macro- 
scopic degeneracy is reminiscent of the degeneracy of the energy levels of a single particle in a 
uniform magnetic field, leading to the well known Landau level structure . An equivalent 
degeneracy occurs, still at the single particle level, for an isotropic 2D harmonic oscillator of 
frequency ujj_, considered in a frame rotating at frequency fl — u± [32| . However we empha- 
size that our result here holds not for a single particle spectrum, but for the eigenmodes of a 
N— body system treated in the mean-field approximation. The occurrence of such a macro- 
scopic degeneracy raises interesting questions concerning the linear response of the rotating 
system to an arbitrary excitation. 

When the quartic potential is present, the above degeneracy is lifted. We set x = r/R in 
with R — (4/i/fc) 1 / 4 , to get the dimensionless eigenvalue equation for the m = modes: 

1 d ( 9 ^ d[<5/°]\ . / x r- i , mR 2 



x (l-eJ-x*)=!p-)=A(e)[6p] where e = ^(cul-n 2 ) . (6) 
x ax \ ax J 2p 

For each f2, i.e. for each e, we are interested in the lowest eigenvalue Ao(e) of the hermitian 
operator in the left hand side of © . We can then deduce the frequency of the lowest transverse 
monopole mode which is plotted as a continuous line in figure |21 



< P ^ + ^£). (7) 



The quantity Aq(c) is an increasing function of e which we calculate numerically using a 
variational method with polynomial trial functions. Just at the critical rotation Q = u>x, 
e = and we get A(0) ~ 11.5. The slow rotation limit corresponds to e ^> 1, in which case 
A(e) ~ 8e; we then recover the result w mp = 2u) j_. 

This analysis also explains the periodic apparition of a hole at the center of the density 
profile for £1 ~ uj± , as observed in figure The stationary density profile ® varies as r 4 
around the origin, whereas the mode Sp varies as r 2 in this region. The curvature of the 
density profile is therefore dominated by the phase of the perturbation. 

For > lu± and a relatively strong excitation (I' /I = 4), we observe that the time 
evolution of the mode structure becomes asymmetric (sec figure |3J this asymmetry is also 
slightly visible in the last row of figure It consists of a periodic entering positive density 
wave, which starts on the edge of the condensate and gradually moves to the center in a time 
period ~ tt/uj±. A possible interpretation of this unusual structure is that two (or more) 
transverse modes m = arc simultaneously excited with a non-zero relative phase. Since 
these modes have similar frequencies (cf. Eq. (JHJ) the initial phase difference between them 
stays nearly constant, and this can give rise to the observed phenomenon. 

To summarize we have presented in this paper a detailed study of the transverse monopole 
mode of a fast rotating degenerate Bose gas. We have shown that the non-harmonic character 
of the potential (which is essential for the confinement of the gas) has a clear influence on 
the mode frequency. We have also given a simple analysis of this mode for an infinitely long 
condensate which is in good agreement with the experimental data. 

Theoretical studies have shown that the addition of a quartic term in the harmonic po- 
tential may lead to the formation of a "giant vortex" , i.e. a vortex with a circulation larger 
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Fig. 3 - Time evolution of the density profile of the oscillating cloud, for f2/27r = 68 Hz. 



than the single quantum h/m .18, 19 20, 21 , 22 33j- The large hole appearing at the center of 
the condensate during the monopole oscillation when f2 > u>± should not be confused with 
such a giant vortex. We are dealing here with a transient state of the condensate, while the 
predicted giant vortex state is a stationary state of the system, for an appropriate angular 
momentum. Another example for a transient large core in a condensate is provided by an 
experiment recently performed in Boulder, in which a hole pierced in a rotating condensate 
confined in a purely quadratic trap was shown to persist for a long time |34| . 

We now briefly discuss some perspectives opened by this work. First we recall that the 
nature of the gas when the rotation frequency Q is around or above u>± is still unclear. The 
detailed study of showed that the number of visible vortices is not sufficient to account for 
the measured rotation frequency of the gas. Two classes of explanation have been proposed to 
account for this finding. Either the fast rotating gas cannot be described anymore by a single 
macroscopic wave function [SSj , or the vortices are still present but they are distorted and do 
not show up clearly in the images of the condensate. We hope that the present experimental 
measurement of the transverse mode frequency can be used to discriminate between these two 
hypotheses. In principle a generalization of the above theoretical treatment to a condensate 
in a 3D trap is possible, and the quantitative predictions of such a model could be compared 
with the experimental findings. 

From the experimental point of view a natural extension of the present work is to switch to 
a two-dimension geometry, using a strong confining potential along the rotation axis. Vortices 
then become point objects and the predicted properties of the system depend on the ratio 
between the atom number N and the number of vortices N v . When N is large, the Bose- 
Einstein condensate presents a regular vortex array, as already observed in ^j- The array 
is expected to melt when N decreases to a value of the order of N v (see also QI]). For 
N < N v the quantum Hall regime for particles should emerge. It is an interesting problem 
to determine the signature of these various regimes on the eigenmodes of the system. For the 
monopole mode considered here and for a pure harmonic confinement, it has been predicted 
in 37 , 38 that the frequency u) mp of the breathing motion of a 2D gas remains strictly equal 
to 2uj± for any equilibrium state. This does not hold anymore in presence of the quartic term, 
and thus the deviation of u) mp from 2 u)j_ can in principle be used to monitor the emergence 
of new quantum phases in the rotating gas. 
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